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INTRODUCTION 

Deardorff's  (1972)  similarity  theory  for  convective  boundary 
layers  (CBL)s  suggests  that  CBL  parameters  are  controlled  by 
buoyancy,  B  =   g/6,  (where  g  is  gravitational  acceleration  and  G 
is  potential  temperature),  surface  temperature  flux,  w'9'0, 
height  z  above  the  surface,  and  zx,    the  height  of  the  inversion 
base  which  marks  the  top  of  the  mixing-layer  (ML).   The  velocity, 
temperature,  and  humidity  scales:  w»,  u»ML,  9»M\  q»ML,  formed  from 
these  parameters  are  used  in  universal  relationships  for  the 
gradients  and  fluxes  of  mean  and  turbulent  quantities  in  the  CBL. 
The  first  of  these,  the  convective  turbulence  scaling  velocity 


w,  =  (Bz,  W6'0)1/3     ,  (1) 

determines  much  of  CBL  physics,  including  diffusion.   Along  with 
w'0'o,  w*  forms  the  basis  for  the  remaining  parameters  used  to 
scale  the  heat,  moisture,  and  momentum,  as  well  as  passive  scalar 
fluxes.   However,  w»  neglects  mechanical  turbulence  due  to  shear 
and  is  only  strictly  valid  for  free  convection  induced  by  surface 
heating.   Yet,  w»  is  used  often  in  mixed  forced/free  convection 
cases  for  lack  of  alternatives.   Here  we  seek  w»s,  an  amendment 
to  w»,   which  includes  mechanical  turbulence  contributions  under 
mixed  forced/free  convection  conditions.   We  assume  mechanical 
generation  comes  from  three  sources:  the  barotropic,  horizontally 
homogeneous  surface  layer  (SL) ,  barotropic  boundary  layer 
entrainment  zone  (EZ) ,  and  a  baroclinity  (BC)  component.   Beyond 
the  z  and  z,  length  scales,  we  propose  that  mixed  forced/free 
convection  also  depends  on  the  SL  roughness  length,  z0,  SL 
Obukhov  length,  L  =  -u»39/ (gkw '9 '0)  (where  u»  is  friction  velocity 
and  k  =  0.4  is  the  von  Karmen  constant),   EZ  thickness,  Ah,  and 
zb,  a  baroclinic  length  scale  which  we  define  below.   We  model 
these  mechanical  sources,  using  relevant  ratios  of  the  length 
scales,  thereby  fulfilling  similarity  requirements. 


APPROACH 

We  outline  the  basic  approach  before  discussing  details.   For 
non-baroclinic,  horizontally  homogeneous,  steady  state  cases,  the 
vertical  profiles  of  the  buoyancy  and  momentum  fluxes  must  be 
linear,  if  the  shapes  of  the  vertical  profiles  of  potential 
temperature  and  wind  speed  are  to  remain  static.   That  is, 


6w»9'z  =  Bw'9'0(l  -   (1  +  a,)z/z,)   ,  (2a) 


-u,'w'z  =  -u,'w'0(l  -  (1  +  a:)z/z,)   ,  (2b) 
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a,  =  -w'e'^/w'G'o  and  a2  =  u,  'w'^/Uj'w'o  (see  fig.  1).   The 
ripts  "zi"  and  "0"  refer  to  conditions  at  the  inversion  base 


urf ace . 


-ii » 


u '  w 


u,' 


w'6'  -   We^e,  and  -U;'w'  •  =  AU,We, 


the  CBL  entrainment  rate,  We  =  zit  and  a6  is  the  potential 
rature  jump  across  the  EZ.  The  subscript  ",t"  denotes  the 
al  derivative  with  respect  to  time.   Such  notation  will 

to  spatial  and  temporal  derivatives  in  general.   The  Betts 
)  correction,  w'e'zi  —    (a6  -  ruAh)We,  hardly  affects  egn.  2a 

Kamada  (1988a)  showed  that  the  much  smaller  interfacial 
ness,  lm,  must  replace  the  EZ  depth,  Ah.   For  lack  of  better 
mation,  we  assume  that  above  zt   the  barotropic  flux  profiles 
inear  with  height  up  to  H  =  zt   +  Ah,  the  top  of  the  boundary 
I.e.. 


6w'6'z  -  Bw'e'zi(l  - 


z  -  z, 


Ah 


(3a) 


■u 


w'  =  -u 

z        I 


w'2i(l  - 


Z   -   Z; 


Ah 


(3b) 


Using  eqns.  (1),  (2a),  and  (3a),  the  vertically  integrated 
buoyancy  flux  or  total  CBL  buoyancy  production  rate  of  turbulence 
kinetic  energy  (TKE)  can  be  written  to  show  that  the  purely  free 
convective  turbulent  scaling  velocity  is  proportional  to  the  cube 
root  of  the  boundary  layer  integrated  TKE  production  rate, 


1  -  a, (1  +  a) 


\       Bw'6^  dz    =   n  w»: 
z0 


where  n  = 


(4) 


and  a  =  Ah/z,.   n  =  0.4  for  mature  CBLs  because  under  such 
conditions,  a,  -  0.2  and  a  ~  0.1.   Since  this  may  not  be  true  for 
other  conditions,  we  discuss  n  in  detail  below.   However,  here  we 
propose  that  this  sort  of  proportionality  extends  to  shear;  that 
mechanical  contributions  to  a  shear  inclusive  turbulent  scaling 
velocity,  w*s,  may  be  written  in  terms  of  the  boundary  layer 
integrated  mechanical  generation  rate  of  TKE.   I.e., 


j   Bw'e'z  dz      +      j-u/w'u.^dz 
z0  zO 


n 


1/3 


(5) 


fig.l    Schematic    profiles    of    a)    buoyancy,    b)    stress,    c)    U    velocity 
and    d)    TKE    shear    generation    rate    for    barotropic    case 
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We  also  propose  that  the  CBL  shear  generation  rate  can  be 
represented  by  discrete  barotropic  surface  layer  (SL) , 
entrainment  zone  (EZ) ,  and  baroclinic  (BC)  components. 


SL                  EZ  BC 

-L H  H 


J-Uj'w'Uj^dz  =   j-u1'w*ulz3z  +   j-Uj'w'u^dz  +   f  -u, '  w  'uib  zdz 

zO  zO  zi  0 

(6) 

Here  uib  z  is  just  the  baroclinic  portion  of  the  vertical  shear. 


BAROTROPIC  SURFACE  LAYER  SHEAR 

We  first  consider  the  barotropic  SL  and  assume  that  most  of  the 
shear  occurs  below  the  height,  -L.   From  SL  similarity  we  have 

u7  =  (u,/kz)0m    ,   with    0m  =  (1  -  15Z/L)'1''3    ,     (7a, b) 


(Carl  et  al.  1973).   Some  workers  prefer  <pm   formulas  involving 
the  -1/4  power.   However,  the  choice  is  among  empirical 
expressions  and  integrals  of  -1/4  power  expressions  result  in 
awkward  complex  numbers.   The  accepted  value  in  the  free 
convection  limit  is  -1/3,  so  we  choose  this  for  theoretical 
appeal  as  well  as  mathematical  convenience. 

Since  0m  is  only  valid  for  z  >>  z0,   we  set  the  upper  integration 
limit  at  -L,  the  lower  limit  at  a  somewhat  arbitrary  10  z0,  and 
treat  the  layer  below  10  z0  separately.   Above  10  z0,  from  eqns. 
(2b)  and  (7a),  we  have 

-l u*3    -l      <Pmdz  u*(u*2   -    AU,We)     -L 

J     -Vw'uiiZdz     =     J        +     1      cpjz 

lOzO  k    lOzO  Z  kz;         10z0 


I 


U,3   -L   0mdz  -L 

-   1    -  (l-a2)  (L/zi2)w,3  {  0mdz 

k   10z0  Z  lOzO     .      (8) 

From  the  L  definition,  u»3/k  =  (-L/z^w,3.   With  the  axes  aligned 
with  the  mean  wind,  we  have  -u'1w'uz  —  -u'w'u^,  since  (vz,  wz)  << 
uz  in  the  mixed  layer.   If  shear  supplies  sufficient  TKE,  10  z0 
<<<  -L.   So  in  eqn.  (8)  term  (2)  -   0.54  (1  -  a2)  (L/Z;)2w»3. 
Letting  p  =  -15z/L,  the  first  integral  reduces  to 


=  V3    ATAN 


V/3(2(p+l)  1/3    +    1) 


m 


1/3 


Ln[  (p+l)-/J+  (p+1)  1M+  1] 


1/3 


Ln[(p+l),,J   -    1] 


-L 
lOzO 


The  upper  limit  evaluates  to  =  1.51.   The  first  two  terms  in  the 
lower  limit  -   1.26,  if  L  >>  150z0.   So  in  eqn.  (8)  the  first 
integral  =  0.25  -  Ln  R,   where  R  =    (1  -  150z0/L)  1/3  -  1  <<  1.   A 
5cm  z0  and  150m  L  would  give  150  z0  =  -0.05L.   Above  10  z0,  this 
leads  to  =  4.36u»3/k  from  term  (1)  in  eqn.  (8).   Oddly,  R  varies 
inversely  with  roughness  length.   But  from  u*  =  ku/(ln(z/z0  -  <pm) 
and  the  definition  of  L,  we  see  that 


■  2,, 3 


nz0/L  =  (ngw,e,/k-uJe)z0(ln(z/z0)  -  <pj 


(9) 


So  for  given  levels  of  u,  w'8',  6,  and  z,  the  ratio  nz0/L  and 
thus  term  (1)  may  vary  much  less  than  z0  itself. 

Below  10  z0  we  neglect  non-neutral  contributions  to  u»  and  assume 
a  constant  value,  -ui'w'uz  =  u»3/10kz0.   Then 


10z0 


j     -u1'w'uz3z  ~  u»3/k 

0 


(10; 


Together,  eqns.  (8)  and  (10)  give  an  estimate  of  the  vertically 
integrated  SL  shear  production  rate, 


-L 


j  -u.'w'u^z  -  (1.25  -  Ln  R)u,3/k   , 

0 

=  (-w»3L/z,)  (1.25  -  Ln  R) 


(lla,b) 


Note  that  we  have  dropped  term  (2)  from  eqn.  (8) ,  since  1)  the 
coefficient  is  an  order  of  magnitude  smaller  than  for  term  (1), 
2)  term  (2)  is  second  order  in  zlf       3)  usually,  -L  <  0.4  zi#  and 
4)  a-,  is  of  order  10"1  -  1;  so  some  further  cancellation  occurs. 

From  eqns.  (2)  and  (3),  if  n  =  0.4,  we  have  the  following  ratio 
of  enhanced  to  unenhanced  TKE  contributions, 


n    -    (1.25    -    Ln   R)L/Z: 


n 


1/3 


12. 


W.  at    (l    +    (3.1    -    2.5Ln    R)  ( -L/  z,)  )  1/Jw» 

(12) 


For  our  previous  example,  w,s  exceeds  w,  by  a  factor  of  =  (1  + 
3  3u»3/w»3) 1/3.   In  this  case,  since  u»/w«  =0.54,  SL  shear 
generation  alone  may  enhance  w»s  by  more  than  80%. 


To  be  useful,  w»s  should  maintain  the  same  proportion  to  the  mean 
turbulent  velocity,  tfui/3,  in  both  the  convective  and  neutral 
limits.   We  test  this  by  noting  for  the  barotropic  case  that  the 
cube  root  of  the  vertically  integrated  production  rate  is  usually 
=  2.4  u,  for  SL  shear.   The  EZ  and  BC  terms  are  not  relevant  for 
the  pure  neutral  case.   From  Grant  (1986)  and  Mason  and  Thompson 
(1987),   we  find  that  aui/3  =  1.3  u»  in  the  layer  averaged  neutral 
boundary  layer.   This  implies  that  aui/ 3  =  0.53  w*s.   In  the  pure 
convective  limit,  the  CBL  averaged  value  for  oul/3    is  ~  0.55  w» 
(Caughey  and  Palmer,  1979  and  Deardorff  et  al.  1980) .   So,  the 
ratio  of  w»s  to  the  layer  averaged  mean  turbulent  velocity  is  the 
same  in  the  convective  and  neutral  limits,  as  required. 


ENTRAINMENT  ZONE  SHEAR 

The  EZ  also  generates  mechanical  turbulence.   To  quantify  it  we 
assume,  moreover,  that  shear  is  constant  across  Ah.  Thus, 


\     -Uj'w'u^dz  =  WeAUf/2 


(13) 


Zl 


Equation  (13)  suggests  the  need  for  a  simple  shear  inclusive 
model  to  estimate  the  entrainment  rate,  We.   For  this  purpose  we 
parameterize  the  TKE  equation, 


storage    buoyancy 


transport 


shear      dissipation 


et  =   6  w'e1  +  (u/e  +  p'e/p) 


u • ' u  ' u  •    -vu   u-  • 


at  the  inversion  base  as, 


(14) 


storage 


buoyancy    transport 


shear 


dissipation 


Cstw»_We/Ah  =s 


-BAGWe      +   Ccw»3/Zi  +   CmAU,2We/Ah    -    Cdw»3/Zj 


(15) 


The  scaling  arguments  are  that  1)  again  w'e'zi  =  -WeAG  and  2) 
large  mixed  layer  eddies  drive  the  TKE  transport  rate  into  the 
EZ.   Thus  TKE  transport  should  vary  with  the  TKE  (oc  w,2)  times 
the  large  eddy  turnover  time  scale,  z,/w».   3)  The  shear 
production  rate  is  modeled  as  the  momentum  flux,  AU,We,  times  the 
local  EZ  shear,  aUj/aIi.   4)  To  maintain  a  steady  state  inertial 
subrange,  the  dissipation  plus  buoyancy  destruction  rates  must 
balance  the  large  eddy  TKE  injection  rate  due  to  transport  and 
shear  production.   So  in  the  neutral  limit,  dissipation  should 
scale  like  transport  with  the  change  rate  of  large  eddy  TKE, 
w»3/Zj.   Since  dissipation  occurs  at  small  scales,  unaffected  by 
buoyancy  anisotropy,  we  retain  the  neutral  form.   For  simplicity 
we  neglect  the  small  TKE  storage  rate  term.   Cm  estimates  vary 
somewhat,  (Pollard  et  al.  1981;  Kato  and  Phillips,  1967,  Tennekes 
and  Driedonks,  1981),  but  we  see  little  reason  for  departures  far 
from  unity,  so  we  assume  Cm  ~  1 . 0 .   We  discuss  Cc  and  Cd  below. 

Tennekes  and  Driedonks  (1981)  presented  a  similar  model,  except 
that  here  we  scale  the  shear  production  locally  with  Ah  instead 
of  z  .   Solving  for  the  entrainment  rate  gives, 


(Cc  -  Cd) 

We  =   w*3/Zi    .  (16) 

3a6  -  AUf/Ah 

Therefore,  from  eqn.  (13)  we  have 

h    (Cc  -  Cd)    Ah 

\  -u,'w'ulzdz  = W,3    ,  (17) 

zi  2(Rb  -  1)    Zj 

where   the   bulk   Richardson   number, 

fi>A9Ah 

Rb   =   —         .  (18) 

au,2 

The  EZ  is  transitionally  turbulent  by  nature.   Thus,  we  expect  Rh 
to  fluctuate  about  some  critical  threshold  turbulence  value,  Rcb . 
Rcb  is  much  larger  than  the  1/4  value  accepted  for  the  stable  SL 
because  EZ  TKE  is  mainly  maintained  by  large  eddy  transport  from 
below.   So  shear  generation  is  only  part  of  the  TKE  injection 
rate,  and  the  dominant  eddies  are  also  much  larger  than  in  the 
transitionally  turbulent  SL.   Thus,  the  EZ's  large  negative 
buoyancy  will  destroy  much  of  the  TKE  during  the  eddy  cascade 
before  it  can  be  dissipated  at  small  scales.   For  these  reasons 


Rcb  may  exceed  unity  in  the  EZ. 

Indeed,  our  simulations  to  steady  state  using  such  an  entrainment 
model  suggest  a  typical  value  of  Rcb  =  1.3.   From  the  results  of 
Turner  (1968),  Deardorff  et  al.  (1980),  and  many  others,  we 
suggest  that  Cc  -  Cd  =  0.25  in  the  convective  limit.   So  in  eqn. 
(17)  the  EZ  shear  contribution  =  0.5  aw,3,  where  a   =   Ah/Zj. 

To  assume  that  a  =  0.1  -  0.2  may  often  suffice.   Then,  ax   —   0.2, 
fl  —  0.4,  and  the  barotropic  EZ  contribution  to  w»s  is  only  a  few 
percent.   But  a   can  exceed  unity  during  rapid  entrainment  into  a 
near-neutral  layer  remaining  from  the  previous  day's  CBL  (Nelson 
et  al.  1989).   During  such  times,  entrainment  models  suggest  that 
aj  may  exceed  40%  (Kamada,  1988a).   Then  n  -   0.1,  and  EZ 
enhancement  of  w»s  may  exceed  0.4.   With  SL  shear  this  may  cause 
w,s  to  double.   In  this  case  a  better  diagnostic  or  prognostic 
estimate  for  a  is  needed.   For  a  diagnostic  estimate,  we  begin  by 
assuming  Rb  =  Rcb  and  solve  eqn.  (17)  for  Ah.   This  leaves  unknown 
aUj  and  a6,  the  inversion  temperature  jump.   a6  is  hard  to 
specify  more  precisely  than  -  +/~  30%  because  z{   and  H  are  not 
that  highly  resolved.   Yet,  we  know  that  a9  grows  with  the 
boundary  layer  growth  rate,  Ht,  into  the  layer  aloft  (with  lapse 
rate  Tu) ,  and  shrinks  with  the  mixed  layer  warming  rate,  9mt. 
So,  using  the  simplest  steady  state  assumption, 


A9t  -  ruHt  -  emt 


=  0 


(19) 


The  mixed  layer  warming  rate  is  given  by  the  surface  heating 
rate,  plus  the  entrainment  rate  of  warm  air  from  aloft,  i.e. 


©m.t  =  (We'0  -  W'8'zi)/Zi     . 


(20) 


So,  if  we  approximate  Ht  as  We(l  +  a),  we  have 


a9  =  w'G'o/We  -  ruz,(l  +  a) 


(21) 


Using  eqns.  (16)  and  (18)  and  combining  constants,  this  results 
in  two  solutions  for  a.   We  show  only  the  physical  root, 


a 


2.3au;2 


ru  r  zc 


1/2 


(22) 


Though  not  explicit,  eqn.  (22)  implicitly  involves  a9  because  it 
presumes  through  eqns.  (18  -  22)  that  a  6  depends  on  ru,  We, 
w'8'o,  and  zx.      This  assumes  that  the  capping  inversion  requires 
finite  potential  temperature  and  windspeed  jumps,  whenever  there 
is  horizontal  wind. 

In  eqn.  (22)  we  must  still  estimate  the  windspeed  jumps  across 
the  EZ.   Wyngaard  (1988)  proposed  the  following  first  order 
estimate.   For  steady  flow  over  a  horizontally  homogeneous 
surface,  the  stress  gradient  is  given  by 


-u'w',  =  f(v  -  vj    ,   and  -v'w'   =  f(uB  -  u)   ,     (23a, b) 


g,  I  iZ      -  V-g 


where  f  refers  to  Coriolis  forcing,  and  the  subscript,  g,  refers 
to  geostropic  winds  (Panofsky  and  Dutton,  1984).  For  barotropic 
flow  with  the  axes  aligned  with  the  mean  mixed-layer  wind, 


-u'w*  z  =  fv  —  (We(u  -  u)  +  u»-)/Zi    ,   and 


-v'w'  z  =  f(u  -  u)  =  We(vg  -  v)/Zj    .  (24a, b) 

The  momentum  balance  becomes 

We(ug  -  u)  +  u*2  =  -fZjV    ,   and 

-We  vg  =  fz,(ug  -  u)    .  (25a, b) 

For  barotropic  cases,  if  au  —  (u  -  u)  and  av  =  (vu  -  v)  across 
the  EZ,  then  the  solutions  to  eqns.  (25a, b)  are  to  first  order, 

au  =  We(u»/fz1)2    ,   and 

av  =  u«2/fz.    .  (26a, b) 


Here,  almost  all  of  the  EZ  shear  results  from  turning  the  wind 
into  the  v  direction.   So  for  a  mature  mid-latitude  BL  with  u,  = 
0.4ms"1,  f  =  lO'V1,  Zj  =  103m,  G  =  300  °C,  and  Tu  =  3  x  10"3  °C  m"1, 
we  have  a  =  0.057.   However,  during  growth  into  the  near  neutral 
remnant  mentioned  earlier,  conditions  may  be  more  like  ru  =  1  x 
10"3  °C  m"1,  and  z,  =  500m.   This  would  result  in  a  =  1.3.   Values 
spanning  such  a  range  were  reported  in  Nelson  et  al.  (1989) . 
They  also  observed  and  prognostically  modeled  the  typical  early 
morning  to  late  afternoon  hysteresis  found  in  plots  of  a   versus 
We/w».   They  noted  that  such  hysteresis  is  absent  from  previous 


diagnostic  expressions  for  the  EZ  depth  which  have  the  form  a  oc 
Ri»"n,  where  Ri»  =   BZ;a6/w«2  and  1/4  <  n  <  1.   No  hysteresis  occurs 
in  such  expressions,  if  we  assume  (see  eqn.  2a)  that  a,,  the 
inversion  to  surface  heat  flux  ratio,  is  constant  because  then 
We/w,  oc  l/Ri,.   Unlike  such  expressions,  eqn.  (22)  shows 
qualitatively  appropriate  hysteresis.   In  fact,  by  combining 
eqns.  (16),  (21),  (22),  and  (26),  we  see  that 


u»4w» 


We/w,  w,9'0Bf2zl3 


(27) 


Far  from  being  constant,  aw»/We  will  typically  increase  initially 
as  surface  heating  and  u»  grow  in  the  young  day-time  boundary 
layer,  but  will  diminish  later  as  z{   becomes  large. 

If  we  require  more  accuracy,  we  can  couple  eqns.  (16,  18,  19,  20, 
and  26)  with 


r 

2u*av  We 

av2a6  t 

1 
1 

f  Zj2a9 

L 

a62 

1 
J 

Ht  =  We  +  (Rcb/S)  | |  (28) 


to  provide  a  complete  entrainment  model.   Here  the  second  term  in 
eqn.  (28)  gives  the  growth  rate  of  Ah  and  thus  a  prognostic 
estimate  for  a. 

More  physics  may  be  added  by  revising  the  We  expression  to 
include  the  small  TKE  storage  rate  term,  Cs[w*2/Ah,  from  eqn.  (16). 
However,  our  present  focus  is  amending  w»,  rather  than  refining 
the  particulars  of  bulk  entrainment  models;  so  we  choose  to  rest 
with  the  above  diagnostic  and  prognostic  formulations  for  a. 

In  either  case,  hysteresis  in  a   versus  We/w,  is  evident,  and  like 
Kamada  (1988b,  figs.  1  and  2),  depends  on  changes  in  a6  and  zx. 
Both  models  include  shear  as  well  as  heating  due  to  entrainment, 
and  avoids  ad  hoc  surface  temperature  distributions,  or  fixed 
values  of  a,.   Also  note  that  eqns.  (16)  and  (21)  lets  us  specify 
a,  which,  together  with  eqn.  (22),  determines  n  from  eqn.  (4). 


BAROCLINIC  SHEAR 

If  the  geostrophic  wind  varies  linearly  with  height  according  to 
the  thermal  wind  relations,  u   =  -(g/fT0)Ty,  and  vg2  =  (g/fT0)Tx, 
then  from  eqn.  (23)  the  baroclinic  stress  profile  must  be 
parabolic,  rather  than  linear  as  in  the  barotropic  case. 
Wyngaard  (1985)  gave  the  following  simple  estimate  for  baroclinic 
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shear  based  on  the  thermal  wind, 

ub.z  =  blUg.z  +  b2Vg,i    f       vb,z  =  b2Ug,z  +  bivgtZ  t  (29a, b) 

where  b,  =  1/  (m2  +  1)  ,  b2  =  m/  (m2  +  1)  ,  m  =  fZj/w*,  and  f  is  the 
Coriolis  forcing.   If  so,  the  stress  can  be  written  as 

b:z  (gz  pfTy  +  2  fT  py) 
2  pfT 


-Uj'w'  =  u*~  +  .  (30) 


So   the   baroclinic    shear   contribution   becomes 

h  gz,(b2Tx   -    bjTy)  (b2gpZj2Ty   +    6pTu,2) 

J  -Ui'w'uiiZdz  =  : 


0  6  pfT"  (31) 

where  terms  in  b,  -  1  =  -m2/ (m2  +  1)  were  considered  neglectably 
small.   The  baroclinic  integral  is  taken  over  the  whole  vertical 
range  0  to  z,  rather  than  -L  to  zs  because  the  SL  shear  analysis 
only  accounts  for  barotropic  effects  based  on  SL  similarity. 
Scaling  analysis  of  egn.  (31)  shows  that  the  squared  and  cubic 
terms  in  zx   are  small.   If  so,  to  first  order,  the  baroclinic 
shear  contribution  from  eqn.  (31)  reduces  to, 


\     -u,'w'uizdz  =  (z,/zb)w,3   ,  (32) 

o 


where  we  can  define  a  baroclinic  length  scale  as, 

fT  w*3  w»3 

zb  =  =   (33) 

u»2g(b2Tx  -  b,Ty)         u,2  ub_z 

Usually,  zb  >>  z,.   But  in  coastal  areas  mesoscale  temperature  and 
pressure  gradients  are  often  large.   For  a  coastal  z,  of  300  m,  a 
3  °C/100  km  horizontal  cross-wind  temperature  gradient,  We  near 
zero,  f  =  10"4  sec"1,  w,  —  0.8  m  s"1,  and  u»  =  0.4  m  s"1,  then  zb  = 
330  m,  and  w»  may  be  enhanced  baroclinically  by  ~  50%. 
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SUMMARY 

In  summary,  we  account  for  surface  layer  (SL) ,  entrainment  zone 
(EZ) ,  and  baroclinic  (BC)  shear.   Each  source  is  given  in  terms 
of  a  dimensionless  ratio  of  the  relevant  length  scales  z0,  L,  Ah, 
Zj,  and  zb.   When  summed  and  multiplied  by  w»,  this  results  in  a 
turbulent  velocity  scale  which  includes  both  buoyancy  and  shear 
generation  of  TKE,  namely, 


fl  -  (1.25  -  Ln  R)L/Zj  +  0.5Ah/z,  +  zjzx 


n 


1  -  a, (1  +  a) 


Here   n  = 


1/3 
w« 


(34) 


a,  =  WeA9/w'e  '0    , 


the  Zj/surface  heat  flux  ratio, 


(Cc    -    Cd) 


We   = 


Ba9     -     AUf/Ah 


w»3/Zj        ,         the   ML   entrainment   rate, 


a9   -  w'9'o/We    -   ruZj(l    +   a)       ,  the   EZ   9    jump, 


a   =   Ah/z,  — 


2.3  (au2    +av2)  1 

+   


ru  b  z,- 


1/2 


,    the    EZ/BL   depth   ratio, 


au  -  We(u»/fz,)-        ,      av  =  u»2/fz;         , 


the    EZ   velocity   jump, 


1/3 


R    s        (1    -    150z0/L)  "J    -    1 


an   SL   shear   generation   term, 


w, 


zb  = 


u»2  u 


b.z 


the  baroclinic  length  scale 


12 


REFERENCES: 

Betts,  A.K. ,  1974:  Reply  to  comment  on  the  paper,  'Non- 
precipitating  cumulus  convection  and  its  parameterization'. 
Quart.  J.  Roy.  Meteor.  Soc,  100,  469-471. 

Businger,  J. A. ,  J.C.  Wyngaard,  Y.  Izumi,  and  E.F.  Bradley,  1971: 
Flux-profile  relationship  in  the  atmospheric  surface  layer.  J. 
Atmos.  Sci.,  28,  181-189. 

Carl,  M.D.,  T.C.  Tarbell,  and  H.A.  Panofsky,  1973:  Profiles  of 
wind  and  temperature  from  towers  over  homogeneous  terrain.  J. 
Atmos.  Sci.  30,  788-794. 

Caughey,  S.J.,  and  S.G.  Palmer,  1979:  Some  aspects  of  turbulence 
structure  through  the  depth  of  the  convective  boundary  layer. 

Deardorff,  J.W. ,  1972:  Numerical  investigation  of  neutral  and 
unstable  planetary  boundary  layers.  J.  Atmos.  Sci.,  29,  91-115. 

Deardorff,  J.W.,  G.E.  Willis,  and  B.H.  Stockton,  1980:  Laboratory 
studies  of  the  entrainment  zone  of  a  convectively  mixed  layer.  J. 
Fluid  Mech.,  100,  41-64. 

Grant,  A.L.M.,  1986:  Observations  of  boundary  layer  structure 
made  during  the  1981  KONTUR  experiment.  Quart.  J.  Roy.  Meteor. 
Soc. ,  112,  825-841. 

Kamada,  R.F.,  1988:  A  fractal  interfacial  entrainment  model  for 
dry  convective  boundary  layers,  part  I:  Model  description.  J. 
Atmos.  Sci.,  45,  17,  2365-2374. 

,  1988:  A  fractal  interfacial  entrainment  model  for  dry 

convective  boundary  layers,  part  II:  Discussion  of  model  behavior 
and  comparison  with  other  models.  J.  Atmos.  Sci.,  45,  17,  2375- 
2383  . 

Kato,  H.  and  O.M.  Phillips,  1969:  On  the  penetration  of  a 
turbulent  layer  into  stratified  fluid.  J.  Fluid.  Mech.  37,  643- 
655. 

Mason,  P.J.  and  D.J.  Thompson,  1987:  Large  eddy  simulations  of 
the  neutral-static-stability  planetary  boundary  layer.  Quart.  J. 
Roy.  Meteor.  Soc,  113,  413-443. 

Mahrt,  L. ,  and  D.H.  Lenschow,  1976:  Penetrative  convection  at  the 
top  of  a  growing  boundary  layer.  J.  Atmos.  Sci.,  33,41-51. 

Nelson,  E. ,  R.  Stull,  and  E.  Eloranta,  1989:  A  prognostic 
relationship  for  entrainment  zone  thickness,  J.  Appl.  Meteorol., 
28,  885-903. 


13 


Panofsky,  H.A.  ,  and  J. A.  Dutton,  1984:  Atmospheric  Turbulence. 
Wiley  and  Sons,  New  York,  p.  133. 

Pollard,  R.T.,  P.B.  Rhines,  and  R.O.R.Y.  Thompson,  1973:  The 
deepening  of  the  wind-mixed  layer',  Geophys .  Fluid  Dyn.  3,  381- 
404. 

Sorbjan,  Z.,  1989:  Structure  of  the  Atmospheric  Boundary  Layer. 
Prentice-Hall,  Englewood  Cliffs,  New  Jersey,  p.  74. 

Tennekes,  H. ,  and  A.G.M.  Driedonks,  1981:  Basic  entrainment 
equations  for  the  atmospheric  boundary  layer.  Bound. -Layer 
Meteorol.,  20,  515-531. 

Turner,  J.S.,  1968:  The  influence  of  molecular  diffusivity  on 
turbulent  entrainment  across  a  density  interface.  J.  Fluid  Mech. 
23,  639-656. 

Wyngaard,  J.C.,  1985:  Structure  of  the  planetary  boundary  layer 
and  implications  for  its  modeling.  J.  Clim.  Appl.  Meteorol.,  24, 
1131-1142. 

,  1988:  Structure  of  the  PBL.  Lectures  on  Air  Pollution 

Modeling,  AMS,  Boston,  1988,  25-29. 


14 


DISTRIBUTION  LIST 


No. copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  VA  22314 

2.  Dudley  Knox  Library  2 
Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

3.  Director  of  Research  (08) 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

4.  Dept.  of  Physics  (PH)  2 
Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

5.  Dr.  Ray  Kamada/061kd  17 
Naval  Postgraduate  School 

Monterey,  CA  93943-5100 

6.  Capt.  David  Struck  (SSD/CLGR)  1 
Los  Angeles  AFB 

Los  Angeles,  CA   90009-2960 

7.  Prof.  Gordon  Schacher  1 
Dean  of  Faculty  and  Graduate  Studies 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

8.  Capt.  Mike  Moss  (AFESC/RDVS)  1 
Tyndall  AFB,  FL   32403-6001 

9.  Lt.  Col.  Stephen  Pryor  (730  SPW)  1 
Vandenberg  AFB,  CA  93437-5000 

10.  Mr.  Steve  Sambol  (730  SPW/DOW)  1 
Vandenberg  AFB,  CA  93437-5000 

11.  Mr.  Randy  Nyman  1 
ACTA 

Vandenberg  AFB,  CA  93437-5000 

12.  Dr.  Don  Ermak  1 
Lawrence  Livermore  National  Lab 

P.O.  Box  808,  L-262 
Livermore,  CA  94551 

13.  Mr.  Bruce  Kunkel  (PL/GPAA)  1 
Hanscom  AFB,  MA   01730-5000 


14.  Mr.  Bill  Boyd  (45  WS/WER) 
Patrick  AFB,  FL   32925-5000 

15.  Mr.  John  Edwards 
SSD/DEV 

Los  Angeles,  AFB,  CA   90009-2960 

16.  Mr.  Larry  Mendenhall 
Geodynamics  Corp. 
21171  Western  Ave 
Suite  100 
Torrence,  CA   90501 

17.  Mr.  Bart  Lundblad  (Aerospace/El  Segundo) 
PO  Box  92957 

Los  Angeles,  CA   90009 

18.  Mr.  Jose  Caraballo 
NASA 

Industrial  Saf ety/RT-SOE 
KSC,  Florida  39899 

19.  Mr.  Jan  Zysko 
NASA 
PM-PCO-4 

KSC,  Florida  39899 

20.  Dr.  Tetsuji  Yamada 
Yamada  Science  and  Art 
147  Monterey  Bay  Dr.  S. 
Los  Alamos,  NM   87544 

21.  Dr.  D.  Pepper 

5301  N.  Commerce  #A 
Moorpark,  CA  93021 

22.  Dr.  Roger  Gibbs 

Naval  Surface  Warfare  Center 
Dahlgren,  Virginia  22448-5000 

23.  Dr.  Kenneth  Moe 
USAF  Space  Division 
Box  92960  World  Way 

Los  Angeles,  CA  90009-2960 

24.  Dr.  Torben  Mikkelsen 
Meteorologie  und  Vind  Energie  Dept. 
RIS0  National  Laboratory 

DK-4000 
Roskilde,  Denmark 


25.  Dr.  Neals  Otto  Jensen 
Meteorologie  und  Vind  Energie  Dept. 
RISO  National  Laboratory 

DK-4000 

Roskilde,  Denmark 

26.  Dr.  John  C.  Wyngaard 
Dept.  of  Meteorology 
Pennsylvania  State  University 
University  Park,  PA  16802 

27.  Dr.  Gary  Geernert 

ONR,  Office  of  the  Chief  of  Naval  Operations 
Arlington,  VA  22217 

28.  Dr.  Ted  Tsui 

Naval  Research  Laboratory 
Monterey,  CA  93943-5006 

29.  Dr.  Steve  Burke 

Naval  Research  Laboratory 
Monterey,  CA  93943-5006 

30.  Dr.  Paul  Gudiksen 

Lawrence  Livermore  National  Lab 
P.O.  Box  808,  L-262 
Livermore,  CA  94551 


DUDLEY  KNOX  LIBRARY 


3  2768  00322869  3 


